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■ The QCD phase diagram is analyzed within chiral quark models with nonlocal 

^ ' separable interactions. For the case of two light flavors, we describe in detail the 

_|_ i characteristics of the chiral phase transition in the T — fi plane, and analyze the 

, color superconducting phase transition taking place at relatively large chemical 

potentials. We also consider the three-flavor version of the models, analyzing the 
predictions for masses and decay constants of light pseudoscalar mesons. 
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The behavior of strongly interacting matter at finite temperature and/or 
density is a subject of fundamental interest, with applications in cosmol- 
ogy, astrophysics, and in the physics of relativistic heavy ion collisions. In 
■ view of the difficulty found when dealing with strong interactions in the 

nonperturbative regime, it has not been possible yet to obtain detailed in- 
formation about the full T — \i phase diagram directly from QCD, and most 
of the theoretical work on the subject relies in low energy effective models. 
One of the most popular ones is that proposed long ago by Nambu and 
Jona-Lasinio (NJL) 1 , in which quarks interact through local four fermion 
couplings. In the last few years, covariant nonlocal extensions of this model 
have been proposed 2 , mainly motivated by some approaches to low energy 
quark dynamics (instanton liquid model, Schwinger-Dyson resummation 
techniques) in which nonlocality arises in a natural way. It has also been 
shown that these nonlocal theories present several advantages in compari- 
son with the standard NJL scheme. Here we analyze, in this context, the 
structure of the QCD phase diagram. First, we consider the case of two- 
flavor models, studying the features of the chiral phase transition for finite 
T and /i and the presence of a possible superconducting phase in the re- 
gion of large /i. Then we also address the three-flavor case, studying the 
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predicted pattern of pseudoscalar meson masses and decay constants. The 
analysis summarized here is developed in detail in Refs. 3,4 ' 5,6 . 

2. Formulation 

We start by quoting the Euclidean quark effective action 



where tp = (u d) T , while the fermion currents j s (x) and j£(x) are given by 



Here r(x — y) is a soft covariant nonlocal regulator, which in the case of 
the instanton liquid model has a definite form in terms of modified Bessel 
functions. In the following we will consider the relatively simple case of 
a Gaussian regulator r(p 2 ) = cxp(— p 2 /2A 2 ), where A is a free parameter 
of the model, playing the role of an ultraviolet cutoff scale. In any case, 
our results are found to be qualitatively similar for other regulator shapes. 
Notice that the constituent quark mass m c (we assume isospin symmetry) 
and the coupling constant G are also free parameters of the model. 

We perform now a bosonization of the partition function of the theory, 
Z = J DtpDip exp(— Se)- To do this we introduce scalar and pseudoscalar 
meson fields a(x), TTk(x), integrating out the quark degrees of freedom. 
Then we carry out the mean field approximation (MFA), in which the meson 
fields are expanded around their translational invariant vacuum expectation 
values a and 7Tfc (the latter vanish due to charge and isospin symmetries). 
The extension to finite T and [i is obtained by performing the replacement 



in the four-momentum integrals arising from the fermion determinant (here 
uj n are Matsubara frequencies corresponding to fermion modes). The grand 
canonical thermodinamical potential per unit volume is then given by 



p 

where S(p 2 ) = m c + a r 2 (p 2 ) is the constituent (momentum-dependent) 
quark mass. Now a(T, fi,m c ) can be obtained from d(jj MFA /da = (gap 
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equation) , and one can easily derive the expressions for other relevant phys- 
ical quantities, such as the quark condensate (qq), quark density p, chiral 
susceptibility Xv,t,h = — (d(qq) /dm c )T,^, pressure, energy density, etc. 

3. Chiral phase transition: analytical and numerical analysis 

In order to study the physical quantities, the Matsubara sums can be 
worked out using complex plane analysis. While in the standard (local) 
NJL model the Euclidean quark propagator has just two simple poles at 
P4 = ±iy / p 2 + M 2 , being M the quark constituent mass, in the case of 
nonlocal models the regulator leads in general to a complicated structure 
of poles and cuts in the complex plane. Depending on the parameters of 
the model and the shape of the regulator, there may be some poles on 
the real axis or not; it has been suggested that this last situation is a sort 
of realization of confinement, in the sense that the effective quark propa- 
gator has no poles at real energies 7 , and one can talk about "confining" 
and "nonconfining" parameter sets. From our analysis it is seen 3 ' 4 that 
using adequate integration paths, and introducing some auxiliary complex 
occupation number functions, the Matsubara sums can be converted into 
integrals. Moreover, we have found that in all considered cases the behavior 
of the physical quantities is strongly dominated by the contribution of the 
first pole of the quark propagator. 

We have carried out numerical calculations to evaluate the dependence 
of the mentioned physical quantities with T and p. To proceed, we have 
adjusted the free parameters A, m c and G so as to reproduce the empirical 
values of the pion mass and decay constant, whereas as the third input we 
have chosen some fixed values of S(0). Both confining and nonconfining 
parameter sets have been considered, studying in each case the features 
of the chiral phase transition. In general, the results have been found to 
be qualitatively similar (a detailed comparison can be found in Refs. 3 ' 4 ); 
for definiteness we present here the curves corresponding to one of these 
parameter sets, namely A = 760 MeV, m c = 7.7 McV and G = 30 GcV~ 2 . 

Our main results are shown in Fig. 1, where we quote the values of the 
quark condensate, the density (relative to nuclear matter density po) and 
the chiral susceptibility as functions of T and p. By looking at the curves 
for the quark condensate, it is seen that for T = there is a first order 
phase transition at a given critical chemical potential p c . The latter goes 
down as T is increased, until one reaches an end point where the first order 
phase transition turns into a smooth crossover. For this parameter set the 
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Figure 1. Behavior of the chiral condensate (left), the density (center) and the chiral 
susceptibility (right) with the temperature and the chemical potential. 

end point is located at (Te,/j,e) = (55 MeV, 210 MeV). In the crossover 
region, the transition is also well defined by the presence of marked peaks 
in the chiral susceptibility, the sharpness of these peaks giving a measure 
of the crossover steepness. Thus we arrive at the phase diagrams shown in 
Fig. 2, in which we plot the phase transition temperature as a function of 
the chemical potential (left) and the density (right). The latter includes a 
"mixed phase" zone, which can be interpreted as a region where droplets 
containing light quarks of mass m c coexist with a gas of constituent, massive 
quarks. It should be noticed that the model predicts a critical temperature 
at yu = of about 100 MeV, somewhat below the results obtained from 
lattice calculations which suggest this value to be ~ 140 — 190 MeV. 




Figure 2. Phase transition temperature as a function of the chemical potential (left) 
and the density (right). Solid (dashed) lines stand for first order transition (crossover). 



4. Extensions to diquark couplings and three flavors 

The ability of this kind of models to allow for a two-flavor superconducting 
(2SC) phase at large chemical potentials is also under study 5 . This analysis 



February 2, 2008 0:50 Proceedings Trim Size: 9in x 6in 



gdummhep 



5 

considers an effective quark-quark interaction in the Euclidean action (1), 
carrying a new coupling constant h. Preliminary results show that for 
low T a first order phase transition into a 2SC phase is expected at some 
critical value of \x. A similar behavior is found when T is increased, until 
one reaches a "triple point" in which all three phases (hadronic, chiral and 
2SC) coexist. The phase diagram is then completed with the inclusion of a 
second order phase transition line between the chiral and the 2SC phases. 
Unfortunately, the precise features of the phase diagram depend on the 
value of the diquark coupling h, which cannot be fixed from low energy 
phenomenology. A detailed discussion in this sense is given in Ref. 5 . 

Finally, we have considered 6 the inclusion of strange degrees of freedom, 
extending the flavor symmetry of the model to SU (3) and incorporating the 
the U(1)a symmetry breaking through a nonlocal 't Hooft-like coupling. As 
a first step we have concentrated on the evaluation of meson masses and 
decay constants, finding our results in quite good agreement with empirical 
values, in particular in the case of the ratio fx/fir and the anomalous decay 
7T° — > 77 (which are problematic in the standard NJL model). We have also 
obtained a reasonable description of the 77—77' phenomenology. Remarkably, 
this requires the presence of two significantly different state mixing angles. 

Clearly, the two-flavor diquark channel analysis cannot be trusted be- 
yond the strange quark mass threshold, and the SU(3) version of the model 
should be addressed, leading presumably to the presence of a color flavor 
locking phase. We expect to report on these issues in forthcoming works. 
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